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Abstract 

A phase diagram for a 2D metal with variable carrier density has been de- 
rived. It consists of a normal phase, where the order parameter is absent; a 
so-called "abnormal normal" phase where this parameter is also absent but 
the mean number of composite bosons (bound pairs) exceeds the mean num- 
ber of free fermions; a pseudogap phase where the absolute value of the order 
parameter gradually increases but its phase is a random value, and finally a 
superconducting (here Berezinskii-Kosterlitz-Thouless) phase. The charac- 
teristic transition temperatures between these phases are found. The chemical 
potential and paramagnetic susceptibility behavior as functions of the fermion 
density and the temperature are also studied. An attempt is made to qualita- 
tively compare the resulting phase diagram with the features of underdoped 
high- Tc superconducting compounds above their critical temperature. 
PACS: 74.72.-h, 74.20.Fg, 74.20.Mn, 64.60.Cn 

I. INTRODUCTION 

The study of the crossover region between superconductivity of Cooper pairs and su- 
perfluidity of composite bosons is attracting much attention due to its close relationship 
to the problem of describing high-temperature superconductors (HTSC) (see, e.g., Refs. 
^^). At present this region is understood for 3D systems, both at zero and finite tem- 
peratures The crossover in quasi-2D systems has also been studied 0, albeit only 
partially, whereas for 2D systems only the case of T = has been studied thoroughly 
0J^. This is related to the fact that fiuctuations of the charged complex order parameter 
in 2D systems are so large that they destroy long-range order at any finite temperature 
(Coleman-Mermin-Wagner-Hohenberg (CMWH) theorem [Q). In this case the appearance 
of an inhomogeneous condensate with a power-law decay for the correlations (the so-called 
Berezinskii-Kosterlitz-Thouless (BKT) phase) is possible. However an adequate mathemat- 
ical description for BKT phase formation is still lacking. 

Most previous analyses [p|-[TT[l of the behavior of 2D systems at T 7^ have been based 



on the Nozieres-Schmitt-Rink approach |12]. This approach is simply a Gaussian approx- 
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imation to the functional integral, and this perhaps explains the difficulties faced in these 
calculations. On the one hand, Gaussian fluctuations destroy long-range order in 2D, and 
if one searches for the T^^ at which order sets in, one should obtain zero in accordance 
with the aforementioned theorems On the other hand, taking Gaussian fluctuations into 
account is completely inadequate to describe the BKT transition |T^. 

Nonetheless, there has been some progress. For example, the BKT transition has been 
studied in relativistic 2 + 1-theory , and the crossover from superconductivity to super- 
fluidity has been considered |jl5[ as a function of the carrier density Uf (see also Ref. [0). 
However, the method employed in Ref. |]T3| to obtain the temperature Tbkt has several 
drawbacks. Most importantly, the equation for Tbkt was obtained without considering the 
existence of a neutral (real) order parameter p, whose appearance at finite T does not violate 
the CMWH theorem. 

As we show below, p defines the modulus of a multivalued complex order parameter $ 
for a 2D system. As a result of allowing for a neutral order parameter, a region where p 
decays gradually to zero appears in the phase diagram of the system. This region separates 
the standard normal phase with p = from the BKT phase, where the correlations exhibit 
power-law decay. Despite the exponential decay of correlations, this new region of states 
may be expected to possess unusual properties, since p plays the same role as the energy 
gap A in the theory of ordinary superconductors in many cases.Q The possible existence of 
such a phase, which in some sense is also normal, may shed light on the anomalous behavior 
of the normal state of HTSC (see, for example, the reviews in Refs. and In 



particular, the temperature dependencies of the spin susceptibility, resistivity, specific heat, 
photoemission spectra, and other quantities [0,|19[ can be explained by the formation of 
either a pseudogap or a spin gap in the region T > T^. 

Using a very simple continuum 2D model, this approach was first attempted in a brief 
note pO|, where we calculated Tbkt and Tp (Tp is the temperature defined by the condition 
p = 0) self-consistently as functions of Uf, and established the boundaries of this new 
pseudogap region, which lies between Tbkt and Tp. 

The purpose of this article is to develop this approach further. Using the static param- 
agnetic susceptibility as an example, we demonstrate that the pseudogap opens below Tp. 
Furthermore, we analyze the difference between the commonly used (see Refs. and [Q) 
pairing temperature Tp and the temperature Tp introduced here. These temperatures turn 
out to be different if the chemical potential p < 0. We also introduce here an abnormal nor- 
mal phase, which lies between Tp and Tp, where preformed bosons exist. This more detailed 
study helps to clarify the physical import of Tp, as well as the nature of the transition at 
Tp. It was believed in the related model that this is a second-order phase transition. 
We argue however, that fluctuations in the phase of the order parameter can transform the 
transition to a crossover, as observed experimentally. 

In Sec. II we present the model and the relevant formalism. The equations for Tbkt? P, 
Tp, and the chemical potentials /x(Tbkt) and p{Tp) are derived in Sec. III. Since the technique 



^ To calculate the observed single-particle spectrum, of course, carrier losses due to scattering of 
carriers by fluctuations of the phase of the order parameter (and in real systems by dopants) must 
be taken into account; see Ref. f[% . 
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employed to obtain the equation for Tbkt is not widely used, we consider it useful to present 
a detailed derivation of this equation. (The details of the calculation of the effective potential 
and useful series are given in Appendix A.) The systems of equations for Tbkt, p(Tbkt), 
/^(Tbkt) and Tp, p(Tp) are analyzed in Sec. IV. The difference between pairing temperature 
Tp and the temperature Tp is discussed in Sec. V. Also discussed is the physical import of Tp. 
Using the example of the static spin susceptibility, it is shown in Sec. VI that the resulting 
pseudogap phase can in fact be used to explain the aforementioned anomalous properties of 
HTSC. 



II. THEORETICAL FRAMEWORK 



The simplest model Hamiltonian density for fermions confined to a 2D volume v is mm\ 



(2.1) 



where x = r, r; ipa{x) is a fermion field, m is the effective fermion mass, /i is the chemical 
potential, and V is an effective local attraction constant; we take Ti = ks = ^■ 



The Hubbard-Stratonovich method, which is standard for these problems can be 
applied to write the partition function Z{v,fi,T) as a functional integral over Fermi fields 
(Nambu spinors) and the auxiliary field $ = l^?/'j?/'J. In contrast to the usual method 
for calculating Z in $, $* variables, the parametrization $(a;) = p{x) exp [— 'j6'(a;)] is more 
appropriate for presenting the corresponding integral in two dimensions p2[ (see also Refs. 
PB| and [0). When this replacement by modulus-phase variables is implemented, it is 
evident that one must also replace ipa{x) = Xo-(a;) exp [i^(x)/2]. Physically, this amounts 
to replacing the charged fermion ipcrix) with a neutral fermion Xcr{x) and spinless charged 
boson e*^(^)/^. Note that while one may formally use any self-consistent definition of the new 
variables, the physical condition that the macroscopic variable $(x) be single- valued under 
2iT rotations fixes the parametrization. This was not taken into account in Ref. where 
a different parametrization was used. 

As a result, one obtains 



Z{v,fi,T) = J pVpV9exp[-pn{v,p,T,p{x),de{x))], 



(2.2) 



where 



pn{v,p,T,p{x),d9{x)) 



dr J dr p^x) -TilnG-^ + TilnGo^ (2.3) 



1 

V 



is the one-loop effective action, which depends on the modulus-phase variables. The action 
( p.3| ) can be expressed in terms of the Green function of the initial (charged) fermions, which 
in the new variables has the operator form 



G 



-1 



-^o'r + rs I — + +rip(r,r) 



7-3 



iV^9{T, r) 
4m 



iV9{T, t)V 



2m 



(2.4) 
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The free fermion Green function Gq = G\f^^p^B=o provides a convenient regularization in the 
process of calculation. It is important that neither the smallness nor slowness of the variation 
of the phase of the order parameter is assumed in obtaining expression ( |2.3| ). In other words, 
it is formally exact. 

Since the low-energy dynamics of phases for which p 7^ is governed mainly by long- 
wavelength fluctuations of 6{x), only the lowest-order derivatives of the phase need be re- 
tained in the expansion of Q{v, fi,T, p{x),d9{x)): 



n{v, p, p{x), dd{x)) ~ ^^kin(^^, Ai, T, p, d9{x)) + VLpot{v, p, T, p), 



where 



(2.5) 



n^^{v,p,T,p,de{x)) = tttJ2 -{gj:r 



n=l 



n 



p=const 



and 



p=const 



(2.6) 



(2.7) 



The kinetic f^kin and potential ^pot parts can be expressed in terms of the Green function 
of the neutral fermions, which satisfies the equation 



and the operator 



6m 



iV^e ^ i\/e{T, r)V' 



Am 



2m 



(2.8) 



(2.9) 



The representation ( |2.5| ) enables one to obtain the full set of equations necessary to find 
^BKT, p(7bkt), and p(Tbkt) at given ep (or, for example, p(T) and p(T) at given T and 
ep)- While the equation for Tbkt will be written using the kinetic part ( |2.6D of the effective 
action, the equations for p(Tbkt) and p(Tbkt) (or p(T) and p(T)) can be obtained using 
the mean field potential (|2.7| ). It turns out that at a phase for which p 7^ 0, the mean- field 
approximation for the modulus variable describes the system quite well. This is mainly 
related to the nonperturbative character of the Hubbard-Stratonovich method, i.e., most 
effects carry over for a nonzero value of p. 

It is clear that the CMWH theorem does not preclude nonzero (p) and, as a consequence, 
an energy gap for fermion x, since no continuous symmetry is broken when such a gap 
appears. Despite strong phase fluctuations in the two-dimensional case, the energy gap in 
the spectrum of the neutral fermion x can still persist in the spectrum of the charged fermion 
ip p2[, even well above the critical temperaturcR We believe that the pseudogap widely 



^ We note that the specific heat experiments [Q] demonstrated the loss of entropy that occms at 
temperatures much higher than Tc. This can be considered indicative of a degenerate normal state, 
consistent with the existence of a nonzero order parameter (p). 
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discussed in high- Tc cuprates might be attributable to the energy gap of a neutral fermion 
introduced in the way described above, so that the pseudogap itself can be considered a 
remnant of the superconducting gap. The condensate of neutral fermions has nothing to do 
with the superconducting transition; the latter is only possible when the superfluid density of 
bosons becomes large enough to stiffen the phase 6{x). The temperature Tp at which nonzero 
(p) develops should be identified in this approach with the pseudogap onset temperature 
PJim . The strategy of treating charge and spin degrees of freedom as independent seems to 
be quite useful, and at the same time a very general feature of two-dimensional systems. 



III. DERIVATION OF SELF-CONSISTENT EQUATIONS FOR Trkt, NEUTRAL 
ORDER PARAMETER, AND CHEMICAL POTENTIAL 

If the model under consideration is reduced to some known model describing the BKT 
phase transition, one can easily write the equation for Tbkt? which in the present approach 
can be identified with the superconducting transition temperature Tc. Indeed, in the lowest 
orders the kinetic term ( |2.6| ) coincides with the classical spin XY-model p5| , pB |, which has 
the continuum Hamiltonian 



n = ^Jdrmr)r. (3.1) 

Here J is the some coefficient (in the original classical discrete XY-model it is the stiffness of 
the relatively small spin rotations) and 6 is the angle (phase) of the two-component vector 
in the plane. 

The temperature of the BKT transition is, in fact, known for this model: 

Tbkt = ^J. (3.2) 



Despite the very simple form^of Eq. ( p.2| ), it was derived (see, e.g., Refs. and p6[) using 



the renormalization group technique, which takes into account the non-single-valuedness of 
the phase 6. Thus, fluctuations of the phase are taken into account in a higher approximation 
than Gaussian. The XY-model was assumed to be adequate for a qualitative description of 
the underdoped cuprates (see also Ref. [^]), and the relevance of the BKT transition 
to Bose- and BCS-like superconductors was recently discussed in Ref. . 



To expand Okin up to ~ (V6')^, it is sufficient to restrict ourselves to terms with n = 1, 2 in 
the expansion ( p. 61) . The calculation is similar to that employed in Ref. , where only high 



densities n/ were considered at T = 0. Thus, to obtain the kinetic part, one should directly 
calculate the first two terms of the series (|2.6| ), which can be formally written = TTr(^E) 
and = |TTr(^S^Il). We note that E has the structure S = t^Oi + IO2, where Oi 



and O2 are differential operators (see (|2.9| )). One can see, however, that the part of E 



proportional to the unit matrix / does not contribute to Hence, 



^An exponentially small correction is omitted here. 
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where 



(2vr) 



5m 



iuj + T3^(k) - Tip 



(3.3) 



(3.4) 



is the Green function of neutral fermions in the frequency-momentum representation, with 
^(k) = £(k) - p and e(k) = k^/2m. 

The summation over the Matsubara frequencies Un = Ti{2n + 1)T and integration over 
k in (|3.3| ) can be easily performed using the sum (|A.7| ); one thus obtains 



^kii = ^ / dr J dr npip, T, p) 



'idrO {Vdf 



6m 



where 



nF{p,T,p{p,T)) 



m 
2^ 



/i2 + p2 + ^ + 2T In 



1 + exp 



T 



(3.5) 



(3.6) 



This has the form of a Fermi quasiparticle density (for p = the expression (|3.6| ) is simply 
the density of free fermions). 

For the case T = p3| , p9| , in which real time t replaces imaginary time r, one can argue 
from Galilean invariance that the coefficient of dfO is rigorously related to the coefficient 
of (VOy. It therefore does not appear in f^kin- wish, however, to stress that these 
arguments cannot be used to eliminate the term (V^)^ from when T 7^ 0, so we must 
calculate it explicitly. 

The Oi term in S yields 



^Si(0i)=2y^ ^^y^^(^ 5:^y t^kTr[6;(zcu„,k)r36?(2a;„,k)r3] 



'id^e [vefV 



+ 



3m 



Using ( |A.11| ) to compute the sum over the Matsubara frequencies, we find that 

^2\ 2 



^gi(Oi) = -| [dr J drK{p,T,p) (^td^O + 



4m 



where 



m 



ir(p,T,p(p,T)) = - 1 + -^^ 



p VjF+f^ 
=: tanh ■ 



2T 



(3.7) 



(3.8) 



(3.9) 



Obviously, the Oi term does not affect the coefficient of (V^^)^. Further, it is easy to 

(2) 

make sure that the cross term involving Oi and O2 in fij^.^^ is absent. Finally, calculations 
of the O2 contribution to ^l^^ yieldQ 



^ Derivatives higher than (V0)2 were not computed here. 
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n^liO^) = T j^T j dr-^ Y.^j rfkk2Tr[e;(za;„,k)/e;(za;„,k)/]^|^. (3.10) 



n=— oo 



Thus, summing over the Matsubara frequencies (see Eq. ( |A.12| )), one obtains 

q!SA02) = ~ dr I dr [ dk (Ve?. (3.11) 

'"^ ' ^0 J USn'm^J 2\/e(k)+p2 ^ ^ ^ 

cosh 

As expected, this term vanishes when T — 0, but at finite T it is comparable with 
Combining ( p.5|) , ( p.ll|) , and ( p.8|) we finally obtain 



Z JO 



^''c/r 1 c/r [ni.(/x,T,p)^9,^ + J(/x,r,p)(V^)' + ir(/x,T,p) (9,^)2] , (3.12) 



where 



J{fi,T,p{^,T)) = ^nF{fi,T,p) - ^ f dx ^ j^g^^g^, 

4m 47r J-a 2T „ / ^2 



cosh"' A / + 



4T^ 



characterizes the phase stiffness and governs the spatial variation of the phase 9{y). One 
can see that our value of the phase stiffness J(T = 0) coincides with the nonrenormalized 
stiffness used in Ref. p7| . 

The quantity J(p, T, p) vanishes at p = 0, which means that above Tp the modulus-phase 
variables are meaningless; to study the model in this region one must use the old variables 
$ and $*. Near Tp one can obtain from (|3.13|) in the high-density limit (see below) 
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J(p ^ e^, T Tp, p - 0) = ^l^ep ^ 0.016^6^. (3.14) 

Direct comparison of ( |3.12| ) with the Hamiltonian of the XY-model ( |3.1| ) makes it possible 
to write Eq. (|3.2| ) for Tbkt directly: 

-J(p,TBKT,P(At,^BKT)) = ^BKT. (3.15) 

Although mathematically this reduces to a well-known problem, the analogy is incomplete. 
Indeed, in the standard XY-model (as well as the nonlinear cx-model) the vector (spin) 
subject to ordering is assumed to be a unit vector with no dependence on T.0 In our case 



^ There is no doubt that in certain situations (for example, very high T) it also can become a 
thermodynamic variable, i.e., one dependent on T, as happens in problems of phase transitions 
between ordered (magnetic) and disordered (paramagnetic) phases when the spin itself vanishes. 
Specifically, for quasi-2D spin systems it is obvious that as one proceeds from high-T regions, a 
spin modulus first forms in 2D clusters of finite size and only then does global 3D ordering occur. 
We note, however, that this dependence was neglected in Ref. |27|, where the nonrenormalized 
phase stiffness J(T = 0) was used to write Eq. (|3.15|) . 
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this is definitely not the case, and a self-consistent calculation of Tbkt as a function of n/ 
requires additional equations for p and /x, which together with (|3.15| ) form a complete set. 

Using the definition ( |2.7| ), one can derive the effective potential f2pot(f, /x, p) (see Ap- 
pendix A). Then the desired missing equations are the condition dQpot{p)/dp = that the 
potential ( |A.1(]| ) be minimized, and the equality v^^dQpot/dp = —nj, which fixes nj. These 
are, respectively 



' ^ ' :tanhV!!^^, (3.16) 



V J (27r)2 2^^2(k)+p2 2T 



nF{p,T,p) = Uf, 



(3.17) 



where np{p^T^p) is defined by 

Equations ( |3.16| ) and (|3.17| ) comprise a self-consistent system for determining the mod- 
ulus p of the order parameter and the chemical potential p in the mean-field approximation 
for fixed T and nj. 

While Eqs. ( p.l6|) and (|3.17| ) seem to yield a reasonable approximation at high densities 
n/, since they include condensed boson pairs in a nonperturbative way via nonzero p, they 
must certainly be corrected in the strong coupling regime (low densities Uf) to take into 
account the contribution of noncondensed bosons (this appears to be important also for Eq. 
( ^.15| ), which determines Tbkt)- The extent to which this alters the present results is not 
completely clear. Previously, the best way to incorporate noncondensed pairs seems to have 
been the self-consistent T-matrix approximation |TI],S(]|-P^, which allows one to account for 
the feedback of pairs on the self-energy of fermions. However, the T-matrix approach, at 
least in its standard form 



10,3 



fails to describe the BKT phase transition, for which 
one must consider the equation for the vertex. On the other hand, in our approach the 
BKT phase transition is realized by the condition (|3.2|) , while an analog of the T-matrix 



approximation in terms of propagators of the p-particle and the neutral fermion % has yet 
to be elaborated. 

The energy of two-particle bound states in a vacuum 



Eb = -2Wexp 



An 



mV 



(3.18) 



(see Refs. and |Q) is more convenient to use than the four-fermion constant V (here 
W is the conduction bandwidth). For example, one can easily take the limits W ^ oo and 
V ^ in Eq. (^.16|), which after this renormalization becomes 



In- 



^/p? + p^ - p 



du ■ 



1 



'M" + — 



exp \ 



fi +1 



(3.19) 



Thus, in practice, we solve Eqs. ( p.l5|) , (|3.17|) , and ( p.l9| ) numerically to study Tbkt as 
function of ?7,/ (or equivalently, of the Fermi energy ep = irnf/m, as it should be for 2D 
metals with the simplest quadratic dispersion law). 
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It is easy to show that at T = 0, the system ( p.l7|) , (|3.19|) transforms into a previously 
studied system (see Ref. ^ and references therein). Its solution is p = \j2\eh\eF and /i = 
— \ei,\/2 + ep. This will be useful in studying the concentration dependencies of 2A/Tbkt 
and 2A/Tp, where A is the zero-temperature gap in the quasiparticle excitation spectrum. 
It should be borne in mind that in the local pair regime (/x < 0), the gap A equals a//?~+~p^ 
rather than p (as in the case /i > 0) 

Setting p = in Eqs. (|3.16|) and (|3.17|) , we obtain (in the same approximation) the 
equations for the critical temperature Tp and the corresponding value of p: 



lnM2 

Tp TT 



tanhu 
du 

u 



(7 = 1.781) 



(3.20) 



T,ln 



1 + exp 



p_ 



(3.21) 



Note that these equations coincide with the system that determines the mean-field temper- 
ature 7<2D)MF ^_ ji^-j jjiiT^^^^'^'^) [|^, evidently as a result of the mean-field approxi- 
mation for the variable p used here. There is, however, an important difference between the 
temperatures T'^^ and Tp. Specifically, if one takes fluctuations into account, T^^^ goes to 
zero, while the value of Tp remains finite. The crucial point is that the perturbation theory 
in the variables p and 9 does not contain any infrared singularities [^,0], in contrast to the 
perturbation theory in $,$*; thus the fluctuations do not reduce Tp to zero. This is why 
the temperature Tp has its own physical meaning: incoherent (local or Cooper) pairs begin 
to form (at least at high enough nf (see Sec. V)) just below Tp. At higher temperatures. 



only these pair fluctuations exist; their influence was studied in Ref. ||35 



IV. NUMERICAL RESULTS 



A numerical investigation of the systems (|3.15| ), ( p.l7|) , (|3.19|) , and (|3.2(]| ), ( |3.21|) yields 
the following results, which are displayed graphically as the phase diagram of the system. 

a) For low carrier densities, the pseudogap phase area (see Fig. 1) is comparable with 
the BKT area. For high carrier densities (ei?>10^|eb|), one easily finds that the pseudogap 
region shrinks asymptotically as 



Tp — Tbkt 
T, 



47; 



(4.1) 



This behavior qualitatively restores the BCS limit observed in overdoped samples. 

b) For ep < (10 — 15)|eb|, the function Tbkt{^f) is linear, as also confirmed by the 
analytic solution of the system (|3.15|) , (|3.17| ), and ( p.l9|) , which yields Tbkt = Re- 
markably, such behavior of Tc{ep) is observed for all families of HTSC cuprates in their 
underdoped region 0,0], though with a smaller coefficient of proportionality (0.01 — 0.1). 
This indicates the importance of including a contribution due to noncondensed pairs in Eq. 
( P.15| ), which defines Tbkt- 
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It has been shown that for optimal doping, the dimensionless ratio eFl\£b\ ~ 3 • 10^ — 10^ 
]. Thus it is quite natural to suppose that in the underdoped region one has e^/lebl ~ 
10 — 10^, where we find linear behavior. 

We note that in this limit, the temperature of formation of a homogeneous order 
parameter for the quasi-2D model |^,^ can easily be written in the form 

T ^ ^^^^^ U 2) 

where t\\ is the interplane hopping (coherent tunneling) constant. This shows that when < 
^BKT, the weak three-dimensionalization can preserve (in any case, at low n/) the regions 
of the pseudogap and BKT phases, which, for example, happens in the relativistic quasi-2D 
model [^. At the same time, as the three-dimensionalization parameter t\\ increases, when 



Tc > Tbkt the BKT phase can vanish, provided, however, that the anomalous phase region 
and both temperatures Tp and Tc ~ n//m are preserved. 

c) Figure 2 shows the values of nj for which /x differs substantially from ep, or in other 
words, the Landau Fermi-liquid theory becomes inapplicable to metals ( also called bad 
ones) with low or intermediate carrier density. As expected, the kink /i at T = Tp, which 



has been observed experimentally and interpreted for the 1-2-3 cuprates |^, becomes 
less and less pronounced a,s ep increases. But in the present case it is interesting that in the 
hydrodynamic approximation employed here, it happens at the normal-pseudogap phase 
boundary or before superconductivity really appears. It would therefore be of great interest 
to perform experiments that might reveal the temperature dependence /i(T), especially for 
strongly anisotropic and relatively weakly doped cuprates. 

d) It follows from curve 3 in Fig. 2 that the crossover (sign change in /i) from local to 
Cooper pairs is possible not only as ei? increases, which is more or less obvious, but also (for 
some Hf) as T increases. 

e) Finally the calculations showed (see Fig. 3) that the ratio 2A/Tbkt is greater than 
4.7 in the region under study. The value 2A/Tp(= 2A/T^^^) is, however, somewhat lower 
and reaches the BCS theory limit of 3.52 only for ep ^ l^fel- It is interesting that this 
concentration behavior is consistent with numerous measurements of this ratio in HTSC 



39| , |40[| . Note that the divergence of 2A/Tbkt and 2A/Tp at ei? ^ is directly related to 



the definition of A at yL^ < 0. 



V. PAIRING TEMPERATURE Tp VERSUS CARRIER DENSITY 

There is no disagreement concerning the asymptotic behavior of Tbkt (or ~ ei? in the 
region of low carrier densities. In contrast, the behavior of the temperature Tp, below which 
pairs are formed, cannot be considered to be generally accepted. For example, in Refs. 
and p^, based on qualitative arguments, this temperature is taken to be the temperature 



Tp of local uncorrelated pairing, which in contrast to Tp increases with decreasing n/.| 



^In fact, in Refs. @ and (see also Ref. [^) this temperature was plotted as an increasing 
function of coupling constant V ^ which for 3D systems corresponds, to some extent, to the carrier 
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Randeria (see Ref. 0] and references therein), to define the pairing temperature Tp, uses 
the system of equations for the mean-field transition temperature and the corresponding 
chemical potential, which is essentially identical to the system ( |3.20| ), ( p.21| ). Thus his 
Tp —i> as n/ 0. 

It is also well known that in the low-density limit, it is vital to include quantum 

fluctuations, at least in the number equation fl^, in the calculation of the critical temper- 
ature at which long-range order forms in 3D. In 2D these fluctuations in fact reduce the 
critical temperature to zero ||Tl|. Certainly quantum fluctuations are also important in the 
calculation of Tp in the limit n j — > and, in particular, in the number equation. However, 
as already stressed in Sec. Ill, these corrections are quite different from what we obtain us- 
ing the variables $, $*, since perturbation theory in the variables p and 6 does not contain 
any infrared singularities p2| , |3^ , and the fluctuations do not yield Tp = 0. In fact, even 
including quantum fluctuations, Tp must exceed Tbkt (p(^bkt) 7^ 0), so that the pseudogap 
phase is always present. 

In our opinion, the temperature Tp has its own physical interpretation: this is the tem- 
perature of a smooth transition to the state in which the neutral order parameter p 7^ 0, and 
below which one can observe pseudogap manifestations. There is also a very interesting and 
important question about the character of the transition. Certainly in the simplest Landau 
theory one appears to have a second-order phase transition, since p takes a nonzero value 



only below Tp 0]. However this kind of transition is only possible for neutral fermions. 
Fluctuations of the ^-phase will transform the pole in the Green function of the neutral 
fermions into a branch cut in the Green function for charged particles in the BKT phase. 
Indeed, the CMWH theorem concerning the absence of spontaneous breaking of a continuous 
symmetry means that symmetry-violating Green functions must vanish. However, it says 
nothing about the gap in the spectrum of excitations, as is sometimes incorrectly stated. 

The correct explanation is that if the symmetry is unbroken, and the fermion excitation 
appears as a pole in the ip two-point function, then the fermion must be gapless. If the 
fermion does not have the same quantum numbers as ip (like our fermion x) and so does not 
appear in the ip two-point function as a one-particle state, then the symmetry does not tell 
whether the fermion (x) will be gapless or not. 

This very general argument ||22| suggests the following plausible scenario. At low temper- 
atures (T < Tbkt), X; P; ^ should be treated as physical quasiparticles (x, p having a gap 
and 6 being a gapless excitation), while a straightforward computation of the two-point 
function reveals its branch-cut structure. 



On the other hand, at temperatures above Tbkt, we should consider ip and $ true 
quasiparticles, since Tbkt is a phase transition point and the spectrum of physical excita- 
tions changes precisely at this point. The ip two-point function at T > Tbkt should be 
studied separately due to the presense of vortices which change the form of the correlator 
{exp[i9{x)]exp[i9{0)]) above Tbkt- In this temperature region the ip two-point function 



loses its branch-cut structure; instead, it acquires the form suggested in Refs. [30| and |31 



density decreasing. In 2D systems, however, where, as is well known, two-particle bound states are 
formed without any threshold, similar conclusions about the behavior of Tp(nj) are questionable, 
and must be checked independently. 
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with a pseudogap originating from the superconducting gap below Tbkt, which preserves 
"BCS-hke" structure as well as the diagonal component of the single-particle Green func- 
tion. In this picture the Fermi-liquid description breaks down, evidently below Tp, due to 
the formation of nonzero p. 

We note, however, that the decisive confirmation of this picture demands further detailed 
study probably based on a different approach, for example the self-consistent T-matrix (see 
Ref. 1^ and references therein), which enables one to directly obtain the full fermion Green 



function. 

To define the temperature Tp properly, one should study the spectrum of bound states 
either by solving the Bethe-Salpeter equation ||^ or by analyzing the corresponding Green 
functions as we do here. It turns out that there is no difference between Tp and Tp in the 
Cooper pair regime (/i > 0), while in the local pair region {p < 0) these temperatures exhibit 
different behavior. 

Indeed, let us study the spectrum of bound states in both the normal (p = 0) and 
pseudogap (p ^ 0) phases. We are especially interested in determining the conditions under 
which real bound states (with zero total momentum K = 0) become unstable. For this 
purpose one can look at the propagator of the p-particle in the pseudogap phase: 



2 5p(r,r)5p(0,0) 



(5.1) 



where pmin is defined by the minimum condition (|3.16|) (or (|3.19|) ) of the potential part (|A.10|) 



of the effective action ( |2.^ ). In the momentum representation, the spectrum of bound states 
is usually determined by the condition 

r^i(a;,K)=0, (5.2) 

where Tp{u), K) is the retarded Green function obtained directly from the temperature Green 
function r(zf2„, K) using the analytic continuation iQ^ —>■ uj + iO. Recall that such analj^ic 
continuation must be performed after evaluating the sum over the Matsubara frequencies. In 
the case of vanishing total momentum K = 0, one arrives at the energy spectrum equation 



From the explicit expression ( p.3|) for Tii{u,0), this function obviously has a branch cut at 
frequencies 

|.| > 2 mmv^5(5)T7 = { ; 1 0° f^'^' 

Thus, bound states can exist below this cut. 

Real bound states decay into two-fermion states when the energy of the former reaches 
the branch point 2min^^^(k) + p^. Since Tp^ is a monotonically decreasing function of tu^, 
it has the unique solution |£6(T)| = 2p(T), at which Eq. ( |5.3|) coincides exactly with the 
mean- field equation ( |3.16D for p(T). It also becomes clear that for /i < we have real bound 
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states with energy eb{T) below the two-particle scattering continuum ai uj = 2a//?~+~P^, 
while at /i > there are no stable bound states. The line /i(T, ep) = in the T-ep plane at 
p 7^ separates the negative /i region where local pairs exist from that in which only Cooper 
pairs exist (positive /i). This line (see Fig. 4) begins at the point T = {e^/7T)\eb\ ~ 0.6|efe|, 
ep ~ 0.39|£b| and ends a.t T = 0, ep = \sb\/2. (The latter follows directly from the solution 
atT = 0,fi = -\et\/2 + ep §0.) 

To find a similar line in the normal phase with p = 0, we consider the corresponding 
equation for the bound states. The propagator of these states (in imaginary time formalism) 
is defined to be 



r-Vr r) = ^'^(^'/^'^.'^(^.r),<|.*(r,r)) 



6^*{t, r)5<l>(0,0) 



(5.5) 

<s>=<s>*=o 



(In the normal phase, where p = 0, we must again use the initial auxiliary fields $ and 
$* (see Sees. II. and III).) Then in the momentum representation (after summing over the 
Matsubara frequencies) we have 



r-'(ffl„,K) = i-i/ 



dk tanh^+(k,K)/2T + tanh^_(k,K)/2T 
e+(k,K)+e_(k,K 




Uk,K) = i-(k±^) -p, (5.6) 



where k is the relative momentum of the pair. The spectrum of bound states is given again 
by Eq. (|5.2| ). Using the energy Sb (see Eq. ( |3.18| )) of the bound state at T = 0, for K = 
we obtain the following equation for the energies of these states in the normal phase: 



dx 



1 tanh(a; -p)/2T' 



\eb\/2 X — jj, — uj /2 



0. (5.7) 



Such states can exist provided — 2p — \eb\ < uj < — 2p. The left-hand side of Eq. (|5.7| ) is 
positive ai UJ = — 2/i — \eb\ and tends to +oo (p > 0) or — oo (p < 0) when uj — 2p. This 
equation always has a solution at p < 0, so bound states with zero total momentum exist 
for negative p. 

For p > 0, analytic analysis becomes more complicated, and requires numerical study. 
One can easily find from (|5.7| ) that at T = 0, stable bound states exist up to p < l^^l/S. In 
fact, numerical study for T > Tp shows that the trajectory fi(T,ep) = (or T = eir/ln2, 
see ( |3.21| )) approximately divides the normal phase into two qualitatively different regions - 
with (p < 0) and without (p > 0) stable (long-lived) pairs. This also holds for other phases, 
which enables one to draw the whole line p(T, ep) = (Fig. 4). 

Knowing the two-particle binding energy, it is natural to define pairing temperature Tp 
as Tp ^ \eb(Tp, fiiTp,ep))\. This equation can be easily analyzed in the region ep <^ \eb\, 
for which we directly obtain Tp ^ \eb\, which clearly coincides with the standard estimate 
PJ^. This means in turn that the curve Tp{ep) starting at Tp(0) ~ \eb\ will be reduced, up 
to the point Tp{0.39ep) ~ 0.6|£b|, which lies on the line Tp(ep) (see Fig. 4). It is important 
that this line is not the phase transition curve; it merely divides the fermion system diagram 
into temperature regions with a prevailing mean number of local pairs (T<Tp) or unbound 
carriers (T>Tp). This is the region of the abnormal normal phase where one has preformed 
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boson pairs. It is widely accepted, however, that this case is only of theoretical interest, 
since there is no Fermi surface (/i < 0) in the phase. The phase area or the difference 
Tp^ep) — Tp{ep) is an increasing function as ep 0, which corresponds to the behavior 



usually assumed pj27 



When /i > there are no stable bound states {eb{T) = 2p(T) = 0) for the normal phase, 
where they are short-lived. Formally, using p(T) = in Eq. ( ^.3|) , we immediately obtain 
( ^.2(j| ) or, in other words, here Tp = Tp. Such a conclusion is in accordance with the generally 



accepted definition of Tp in the BCS case |41 



Thus the phase diagram of a 2D metal above T^. acquires the form shown in Fig. 4. It 
is interesting that if the line Tp{ep) cannot be defined exactly, the temperature Tp{ep) is 
the line below which pairs reveal some signs of collective behavior. Moreover, aX T < Tp 
one can speak of a real pseudogap in the one-particle spectrum, while in the region Tp < 
T < Tp only strongly developed pair fluctuations (some number of pairs) exist, though they 
probably suffice to reduce the spectral quasi-particle weight, and to produce other observed 
manifestations that mask pseudogap (spin gap; see Ref. |^) formation. 



VI. PARAMAGNETIC SUSCEPTIBILITY OF THE SYSTEM 



It would be very interesting to study how a nonzero value of the neutral order parameter 
affects the observable properties of the 2D system. Does this really resemble the gap opening 
in the traditional superconductors, except that it happens in the normal phase? Or, in other 
words, does the pseudogap open? 

We shall demonstrate this phenomenon, taking the paramagnetic susceptibility of the 
system as the simplest case in point. To study the system in the magnetic field H one must 
add the paramagnetic term 



HPM = -^^bH U}(r)^T(r) - i^l{r)^l{r) 



(6.1) 



to the Hamiltonian (|2.1| ) where fxp = eh/2mc is the Bohr magneton. Note that, using the 
isotropy in the problem, we chose the direction of field H to be perpendicular to the plane 
containing the vectors r. 

Adding the corresponding term to Eq. (|2.8| ) for the neutral fermion Green function, it 
is easy to show that in the momentum representation (compare with 



g{iuJn,k,H) 



[iUn + ^^BH)i + r3^(k) - Tip 

(zcu„ + /iBi/)2-e(k)-p2 ■ 



(6.2) 



The static paramagnetic susceptibility can be expressed in terms of the magnetization 

dM{p,T,p,H) 



OH 



(6.3) 



H=0 



which in the mean-field approximation can be derived from the effective potential: 

l9fipot(f,p,T,p, iJ) 



M{p,T,p,H) 



OH 



(6.4) 
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Thus from ( |6.4|) one obtains 

°° r dk 

Then using the definition (| 



(6.5) 



xif^,T,p) = Pb 



one arrives at 
dk 



2T y: 



n212 ' 



(6.6) 



The sum in 
final result 



^ can easily be calculated with the help of Eq. ( [A.ll| ); thus, we obtain the 



1 r 

XiP,T,p) = XPauli^ J_ 



dx 



cosh^ 



(6.7) 



lx' + 



4T' 



where XPauU = p'sm/n is the Pauli paramagnetic susceptibility for the 2D system. 

To study X as a function of T and Uf {or ep), Eq. ( |6.7| ) should be used together with 
Eqs. ( FTtD and (|3l9| ). 

For the case of the normal phase (p = 0) one can investigate the system analytically. 
Thus (|6.7|) takes the form 



X{.P^,T,p = 0) = XPauli- 



'1 + exp(-/i/T) ' 
where fi is determined by ( |3.21| ). This system has the solution 

x(eF,T,p = 0) = XPauii[l - expi-ep/T)], 



(6.8) 



(6.9) 



which is identical to a solution known from the literature 42 



The results of a numerical study of the system (|6.7|) , ( |3.19| ), and ( |3.17|) are presented 
in Fig. 5. One can see that the kink in x occurs at T = Tp as in the dependence of /i on 
T. Below Tp the value of x{T) decreases, although the system is still normal. This can be 
interpreted as a spin-gap (or pseudogap) opening. The size of the pseudogap region depends 
strongly on the doping {ep/\eb\), as observed for real HTSC P|,|TB|,|TP| . For small values of 



Sp/lebl this region is very large (Tp > 6Tbkt), while for large ep/\eb\ 
larger than the region corresponding to the BKT phase. 



5 — 30 it is slightly 



VII. CONCLUSION 

To summarize, we have discussed the crossover in the superconducting transition between 
BCS- and Bose-like behavior for the simplest 2D model, with s-wave nonretarded attractive 
interaction. 

While there is still no generally accepted microscopic theory of HTSC compounds and 
their basic features (including the pairing mechanism), it seems that this approach, although 
in a sense phenomenological, is of great interest since it is able to cover the whole range of 
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carrier concentrations (and thus the whole range of couphng constants) and temperatures. 
As we tried to demonstrate, it enables one to propose both a reasonable interpretation for 
the observed phenomena caused by doping and to describe new phenomena — for example, 
pseudogap phase formation as a new thermodynamically equilibrium normal state of low- 
dimensional conducting electronic systems. 

Evidently there are a number of important open questions. They may be divided into 
two classes: the first concerns the problem of a better and more complete treatment of the 
models themselves. The second class relates to the extent to which this model is applicable 
to HTSC compounds, and what the necessary ingredients are for a more realistic description. 

Regarding the microscopic Hamiltonian as a given model, our treatment is obviously 
still incomplete. In particular, there exists an unconfirmed numerical result |^ based on 
a fully self-consistent determination of a phase transition to a superconducting state in a 
conserving approximation, which states that the superconducting transition is neither the 
simple mean-field transition nor the BKT transition. (See, however, the discussion preceding 
Eq. (|3.18| ).) Besides, it would be very interesting to obtain the spectrum of both the 
anomalous normal and pseudogap phases. It is important also to take into consideration the 
effects of noncondensed bosons, which might help to obtain a smaller slope in the dependence 
of Tbkt on ep. 

As for the extent to which the models considered are really applicable to HTSC, most 
of the complexity of these systems is obviously neglected here. For example, we did not 
take into account the indirect nature of attraction between the fermions, d-wave pairing, 
inter-layer tunneling, etc. Nevertheless, one may hope that the present simple model can 
explain the essential features of pseudogap formation. 
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APPENDIX A: CALCULATION OF THE EFFECTIVE POTENTIAL 



Here we sketch the derivation of the effective potential. To obtain it one must write Eq. 
in the momentum representation: 

fipot(t;,/x,T,p) = t;||,-T Y: y ^Tr[ln^-i(^u;„,k)e^^--3] 

+ ^ E y (^Tr[lnGoi(zc.„,k)e^^--^]|, 6-^+0, (A.l) 
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where 



Q ^{iujn, k) = iuj - r3^(k) + np, Gg k) = ^ k) 



p=^l=o 



(A.2) 



are the inverse Green functions. The exponential factor e^^^"'^^ is added to ( [A.l|) to provide 
the correct regularization which is necessary to perform the calculation with the Green 



functions [44|. For instance, one obtains 

+00 



+00 



liin^ E Tr[ln^-i(icu„,k)e^'^'^""^] = ^lim <j ^ Tr[ln^-i(itu„, k)] cos5w„+ 



5^+0 



i ^ sin5a;„Tr[(ln^ k) - In^ ^(-zu;„, k))r3] 

li'n>0 



X: Tr[ln6;-i(za;„,k)]- 



i(k) 



(A.3) 



where 



^ ^(ia;n,k) -r3^(k) + rip 
m ^ ~ : , ^ 00 



and 



E 



sin (JciJ^ 



aj„>0 



1 /""o sin5x 1 ^ 

= — sign 0. 

X AT ^ 



2tiT Jo 



To calculate the sum in ( |A.3D , one must first use the identity Tr In A = Indet A, so that 
( |A.1| ) takes the form 



f 2 +00 
Qpotiv, fi,T, p) =v I — -T J2 J 



dk , det^ ^(ici;„,k) 
In ■ 



(27r)2 detGo^(«o^n,k) J {2n)^ 



dk 



-e(k)+e(k)]|. 

(A.4) 



Calculating the determinants of the Green functions (|A.2|) one obtains 

f^pot T, p) = . I - - T y ^ In -j^^^^ J [-^(k) + em | , 

(A.5) 

where the role of Go{iuJn,k.) in the regularization of Qpot is now evident . The summation 
in ( |A.5|) can be done if one uses the representation 



In^ 



dx 



1 



u;2 + 52 Jq I l<j2 + + X +h'^ + x ) '' 



(A.6) 



and then 
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°° 1 Vr TTC 

-r--, ?r = — tanh — . 

^0 {2k + 1)2 + c2 4c 2 



(A.7) 



We find 
In 



+ 62 



2v^ 



tanh 



X 



2T 



2Va2 + X 



tanh 



-\/a2 + X ' 
2T 



(A.8) 



Integrating (|A.8|) over x, one thus obtains 
Finally, substituting ( |A.9|) into (|A.5| ), 



rfk , cosh[v'e'(k)+p2/2T] 
m 



(27r)2 



cosh[£(k)/2T] 



(A.9) 



V 



(2vr) 



cosh[Je'(k)+p2/2r] 



(A. 10) 



It is easy to show that at T = 0, the expression ([A.10| ) reduces to that obtained in Ref. [0. 

Finally, we give formulas for the summation over the Matsubara frequencies used in 
Sees. Ill and VI: 



T Tr[^(it^„,k)r3^(it^„,k)r3] = 2T ^ 



[e^(k)+p2]3/2 



tanh 



e(k) + p 
2T 



n=—oo 
2 ^2 



e^(k)-p2-^2 

:a;2+e(k)+p2]2 

e^(k) 1 



2T[e2(k)+p2 



cosh 



e ^(k)+p2 
2T 



(A.ll) 



e2(k) + p2-c.2 



T E Tr[^(...„,k)/^(..;„,k)/] = 2T r, ,2 , .2^^ , ,212 

n=-oo n=-oo + ? l^J + P J 

1 1 



2T 



cosh 



2 ve(k)+p2 



(A.12) 



2T 



where the Green function Q{iujn, k) is given by ( |3.4| ). Both formulas can easily be calculated 
using Eq. ( |A.7|) and its derivative with respect to c. 
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FIGURE CAPTION 



FIG. 1. 

^BKT and Tp versus the noninteracting fermion density. Dots represent the function p(ei?) 
at T = Tbkt- The regions of normal phase (NP), pseudogap phase (PP), and BKT phase 
are indicated. 



FIG. 2. 

/x(T) for various values of eF/\eb\: 1) 0.05; 2) 0.2; 3) 0.45; 4) 0.6; 5) 1; 6) 2; 7) 5. (For 
/J, > and II < the chemical potential was scaled to tp and \6b\, respectively.) The thick 
hnes delimit regions of BKT, pseugogap (PP), and normal (NP) phases. 

FIG. 3. 

2A/Tbkt and 2A/Tp versus the non-interacting fermion density. 
FIG. 4. 

Phase diagram of the 2D metal at low concentrations. The dotted line corresponds to /i = 0, 
and the temperature Tp separates abnormal normal phase (ANP) from normal phase. The 
critical temperature Tbkt is not shown. 

FIG. 5. 

X(T) for various values of ep/leb]- 1) 0.6; 2) 1; 3) 5; 4) 10; 5) 30. 
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FIG. 4. 
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